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1.  INTRODUCTION 

The  purpose  of  this  technical  memorandum  is  to  present  the  R/K  optical  profile  function  and 
several  graphical  comparisons  that  illustrate  the  power  and  versatility  of  the  function  to  fit 
nonlinear  curves.  Here,  we  did  not  review  or  present  any  of  the  multitude  of  mathematical 
curve  fitting  methods  that  are  presented  in  numerous  texts  and  papers.  Instead,  we 
focussed  on  the  R/K  functional  form  that  we  found  to  fit  most  of  the  vertical  profiles  of 
extinction  and  backscalter  coefficients  that  were  simulated  by  our  theoretically  based 
microphysics  model  (Rachele  and  Kilmer,  1992).  The  R/K  function  did  not  evolve  without 
considerable  effort.  We  had  tried  fitting  extinction  (o^)  and  backscatter  (fi)  profile  curves  with 
other  standard  and  nonstandard  functional  forms  and  approaches,  but  we  were  not 
consistently  successful.  One  of  these  approaches  involved  calculating  <7g(z  +  h)  versus 
<7g(z)  data,  where  z  is  height  above  ground  level  and  h  is  a  constant  height  increment,  as 
proposed  by  Duncan  et  al.  (1900).  Those  <ig(z  +  h)  versus  were  transformed  (by 

rotation  and/or  translation  of  axes),  and  fitting  functions  were  investigated  for  the 
transformed  data.  Another  approach  involved  two-piece,  double-exponential  functions 
(Heaps,  1982;  Rachele  and  Kilmer,  1991;  Kilmer  and  Rachele,  1992a)  such  as 

Og  =  A  exp  [B  exp  (Cz)]  ,  (1) 

where  A,  B,  and  C  are  constants  with  appropriate  values  for  the  case  being  approximated. 
The  latter  approach  was  one  of  the  more  promising,  but  was  not  acceptable  for  our 
purposes.  This  was  especially  true  for  simulated  backscatter  coefficient  profiles  where  the 
structure  was  too  extensive  to  allow  successful  fitting  by  this  approach.  Specifically,  we 
were  searching  for  a  continuous  function  that  would  result  in  a  good  fit  over  a  wide  range 
of  points,  while  being  constrained  to  fit  the  end  points  precisely.  In  particular,  we  wanted  a 
function  that  required  a  minimum  number  of  fitting  constants,  and  we  did  not  want  to  fit  the 
entire  curve  by  cutting  it  into  several  subsegments.  For  these  reasons,  we  did  not  consider 
fitting  with  splines  to  be  a  desirable  option. 


2.  THE  R/K  OPTICAL  PROFILE  FUNCTION 


The  R/K  optical  profile  function  is  a  higher  order  exponential  written  as 


where 


(2) 


(3) 


In  equations  (2)  and  (3).  x  is  the  independent  variable  and  y  is  the  dependent  variable.  The 
constraint  points  (end  points)  are  designated  as  (x,,  y,)  and  (Xg,  y2)-  N(x),  a  function  of  x,  is 
represented  in  general  form  as 


N(x)  =  exp 


or,  for  curves  that  contain  oscillations, 

Na 


N(x)  =  Aq  -F  2  ^  An  cos 


n  =  1 


27rn(x  — 


Xp2  -  X 


— 1  -if 


Bn  sin 


27rn(x  —  Xpi) 

“  Xpi  -F  1 


(4) 


(5) 


in  the  range  Xp,  <  x  <  Xpj-  In  equations  (4)  and  (5),  A^,.  A,  through  Ana.  through 

Bnb.  Xpi.  Xp2,  ag.  and  a^  through  an  are  constants.  In  our  work  we  have  found  that  n  =  3  is 
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generally  sufficient  when  using  equation  (4)  to  fit  extinction  coefficient  profiles;  whereas, 
when  fitting  backscatter  profiles,  as  many  as  31  ( =  1  +  +  Nb)  constants  are  used  in 

equation  (5). 

To  implement  ttie  functional  structure  defined  by  equations  (2)  through  (5),  we  first  solve 
equations  (2)  and  (3)  for  N(x)  giving 


N(x)  = 


Using  some  sets  of  N{x)  data  calculated  according  to  equation  (6).  In  N(x)  versus  x  data  were 
fitted  with  polynomials  to  determine  values  of  the  constants  in  equation  (4).  Some  other  sets 
of  N(x)  data  calculated  according  to  equation  (6)  were  fitted  using  Fourier  analysis,  resulting 
in  the  fitting  expression  given  by  equation  (5).  For  some  cases,  one  must  use  a  combination 
of  polynomial  and  Fourier  terms.  The  polynomial  and  Fourier  constants  were  determined 
by  using  IMSL  (International  Mathematical  and  Statistical  Library)  routines.  All  of  the  points 
(x,  y)  selected  for  fitting  must  be  such  that  x.,  <  x  <  Xj.  All  points  used  for  polynomial  fitting 
of  In  N(x)  must  have  y  values  such  that  y^  <  y  <  y?. 

3.  EXAMPLES  AND  RESULTS 

Values  of  (x,.  y,)  simulated  by  three  different  physical/mathematical  models  were  considered 
in  this  study.  For  the  first  example,  values  (y,)  of  extinction  and  backscatter  coefficients 
simulated  with  the  R/K  microphysics  model  were  considered.  A  plot  of  each  of  these 
simulated  profiles  is  shown  as  solid  lines  in  figures  1  and  2,  where  z(  =  x)  is  the  height 
above  ground.  For  simulating  data  used  in  figure  1,  model  input  included  reference  height 
(2  m  above  ground  level)  values  of  95  percent,  5  km,  and  0  °C  for  relative  humidify,  visibility, 
and  ambient  air  temperature,  respectively.  Near  adiabatic  conditions  and  a  maritime  air 
mass  were  also  selected.  Model  input  for  simulating  data  used  in  figure  2  also  specified  a 
maritime  air  mass,  but  a  maximum  liquid  water  content  that  would  be  75  percent  of  the  near 
adiabatic  value  was  selected,  and  reference  height  values  of  99  percent,  2  km,  and  30  °C 
were  specified  for  relative  humidity,  visibility,  and  ambient  air  temperature,  respectively. 

These  curves  were  considered  as  two  segments  for  each  case.  The  first  segment 
represents  the  subcloud  region  for  2  m  <  z  <  103  m  for  extinction  and  2  m  <  z  <  36  m  for 
backscatter.  The  second  segment  represents  the  cloud  region  for  104  m  <  z  <  463  m  and 
37  m  ^  z  <  460  m  for  extinction  and  backscatter.  respectively.  The  fitting  curves  produced 
by  equation  (2)  for  extinction  and  backscatter  coefficients  are  shewn  as  dashed  lines  in 
figures  1  and  2.  We  considered  the  fits  in  each  case  to  be  very  good.  Values  of  (the 
coefficient  of  determination)  were  0.99986  and  0.99855  for  the  fits  shown  in  figures  1  and  2, 
respectively.  The  In  Nf^)  for  extinction  was  fit  with  a  cubic  polynomial,  figure  3.  The  fitting 
curve  for  backscatter  shown  in  figure  4  was  calculated  by  fitting  In  N(z)  with  a  cubic 
poiynomial  for  the  first  segment  and  fitting  N(z)  with  a  16-term  Fourier  series  for  the  second 
segment.  Note  that  even  though  the  N(z)  curves  in  figures  3,  4,  and  6  are  discontinuous  at 
the  cloud  base,  the  corresponding  and  /?  curves  are  continuous  as  shown  in  figures  1,  2, 
and  5.  As  expected,  the  and  /?  curves  pass  precisely  through  the  specified  values  at  the 
reference  height  (2  m),  the  cloud  base,  and  the  cloud  top.  Overall,  the  function  fits  of  o*  and 
/?  are  very  good.  We  did  find,  however,  that  for  some  cases  a  large  number  of  Fourier  terms  • 

were  required  to  obtain  a  good  fit.  As  an  example,  31  Fourier  terms  were  used  for  the 
second  segment  shown  in  figures  5  and  6. 


Ill 


In  y  —  In  y. 
In  yj  -  In  y^ 


In 


(  “  -  'I 

1  ~  y-  / 


(6) 
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For  the  second  example,  we  used  data  produced  by  our*  thermal  conductivity  model  for 
loam  soil,  where  l<-s  (  =  y)  is  the  thermal  conductivity  and  </>vv  (  =  x)  is  the  volume  fraction  of 
water.  When  generating  the  solid  line  curve  shown  in  figure  7,  we  assumed  that  the  volume 
frachon  of  solid  material  (sand,  silt,  and  clay)  was  0.5  and  the  volume  fraction  of  organic 
matter  was  0.10.  This  curve  was  fit  by  equation  (2)  for  0.2  <  kj  <  2.1  and  0  <  4)^  <  0.27. 
Casual  inspection  shows  that  the  fit  is  excellent.  The  In  N((/),^)  for  this  case  was  fit  by  a  fourth 
degree  polynomial,  figure  8. 

The  third  example  makes  use  of  a  curve  produced  by  our  atmospheric  planetary  boundary 
layer  model.**  The  solid  line  in  figure  9  shows  the  change  of  wind  direction  a  ( =  y)  as  a 
function  of  height  z  (  =  x),  where  a  varies  from  13°  to  40°  for  2  m  <  z  <  1000  m.  The  R/K  fit 
to  this  curve  is  shown  in  figure  9.  Again  the  fit  is  excellent.  The  In  N(z)  values  for  this  case 
were  fit  by  a  fourth  degree  polynomial,  figure  10. 

4.  CONCLUSIONS 

For  the  examples  considered  in  this  study,  we  found  that  the  R/K  optical  profile  function  fits 
smooth  model  simulated  curves  very  well.  In  a  separate  study  (Kilmer  and  Rachele,  1992b), 
we  fitted  several  hundred  extinction  and  backscatter  coefficient  profiles  with  the  R/K  optical 
profile  function  and  reached  the  same  conclusion.  Some  cases,  however,  did  require  more 
terms  than  we  would  have  liked  when  fitting  N(z),  o^,  and  /?.  Some  backscatter  coefficient 
curves  with  many  large  amplitude  oscillations  are  examples  of  such  cases.  Some  cases 
required  both  polynomial  and  Fourier  terms. 


*  Rachele  and  Manguso,  U  S  Army  Research  Laboratory,  White  Sands  Missile  Range,  New  Mexico.  88002. 
(Unpublished) 

**  Rachele,  Hansen  and  Manguso.  U.S  Army  Research  Laboratory.  White  Sands  Missile  Range,  New  Mexico. 
88002. 
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Figure  1.  Solid  line:  R/K  microphysics  model  simulation  for  a  wavelength  of  4.0  pm; 

dashes:  calculated  using  R/K  optical  profile  function  as  defined  by  equations  (2). 
(3).  and  (4),  with  x  =  height,  y  =  extinction  coefficient,  and  n  =  3. 


Figure  2.  Solid  line:  R/K  microphysics  model  simulation  for  a  wavelength  of  10.6  pm; 

dashes:  calculated  using  R/K  optical  profile  function  as  defined  by  equations  (2), 
(3),  and  (4)  (tor  x  <  36  m)  or  (5)  (for  x  ^  37  m),  with  x  =  height,  y  =  backscatter 
coefficient  in  equation  (2),  =  8.  and  Ng  =  7. 
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In  N(z) 

Figure  3.  Solid  line:  R/K  microphysics  model  simulation  for  a  wavelength  of  4.0  /jm; 

dashes:  calculated  using  equation  (4),  with  x  =  height,  y  =  extinction  coefficient 
in  equation  (2),  and  n  =  3. 


Figure  4.  Solid  line:  R/K  microphysics  model  simulation  for  a  wavelength  of  10.6  ^m; 

dashes:  calculated  using  equation  (4)  for  x  <  36  m  and  equation  (5)  for  x  ^  37 
m,  with  x  =  height,  y  =  backscatter  coefficient  in  equation  (2),  N;^  =  8,  and 
Nb  =  7. 
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Figure  5.  Solid  line:  K  microphysics  model  simulation  for  a  wavelength  of  5.0 

maritime  air  mass,  near  adiabatic  conditions,  and  reference  height  (2  m)  values 
of  95  percent,  5  km,  and  30  °C  for  relative  humidity,  visibility,  and  temperature, 
respectively.  Dashes:  calculated  using  R/K  optical  profile  function  as  defined  by 
equations  (2),  (3),  and  (4)  (for  x  <  129  m)  or  (5)  (for  x  >  130  m),  with  x  =  height, 
y  -  backscatier  coefficient,  and  =  Ng  =  15. 


N(z) 


Figure  6.  Solid  line:  R/K  microphysics  model  simulation  for  same  wavelength  and  case  as 
in  figure  5.  Dashes:  calculated  using  equation  (4)  for  x  ^  129  m  and  equation  (5) 
for  X  ^  130  m,  with  x  =  height,  y  =  backscatter  coefficient  in  equation  (2),  and 
Nrt  =  NB=15. 
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Figure  7.  Thermal  conductivity  vs.  water  content  for  loam  where  <1)^  =  0.5,  (^q  =  0.2 

\\\  fit  with  RK  optical  profile  fit  N(^y^)  =  exp  (0.040407  +  3.3972  -  17.454 

-  769.18  <f>l  +  2460.47  </)^. 
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